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We consider interaction effects related to a nonuniform magnetization in ferromagnetic nanowires and their possible 
generalization to nanorings. First we show that the electric current in a ferromagnetic nanowire with more than one 
domain wall induces an exchange coupling between the walls mediated by the spin-dependent interference of scattered 
carriers. This interaction reveals a complex behavior as a function of mutual orientations and separation of the domain 
walls, thus affecting the domain wall dynamics. Then we consider how the theory should be modified in the case of a 
magnetized quantum ring. In this case the situation is more complicated because the vortex and onion magnetization 
states in the ring are not truly ferromagnetic. 
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1 Introduction Domain walls (DWs) are noncollinear 
magnetization regions separating areas of different ho- 
mogenous magnetizations. In nanostructures they are the 
subject of numerous investigations and are also important 
for applications [1,2.3]. In low dimensional nanostruc- 
tures, such as the wires and rings we are here concerned 
with, the coupling between the electrons and the DWs is 
enhanced in comparison with the bulk scenario [4|. This 
can drastically alter the transport properties of the nanos- 
tructures 1516 )71. The spin flip of the carriers as they tra- 
verse the region of non-collinearity leads to a spin-torque 
acting on the DW, and consequently to the possibility of 
current-induced DW motion 0121 . 



DWs can exist not only in magnetic nanowires but can 
also be created in magnetic nanorings [8][9] between the 
regions of a locally defined homogenous magnetization. 
Fully ferromagnetic rings, where the magnetization direc- 
tion is globally preserved, are not experimentally found, 
though they inform the majority of theoretical considera- 
tions. Instead, the magnetization direction usually follows 



the shape of the annulufl In this case it is clear that the 
DWs always come in pairs and as such a natural question 
to ask is how they interact, mediated by the electrons in 
the system. As a precursor to this question we discuss such 
interactions, and their effect on the DWs in nanowires 1 10 
[TTI . This will form the bulk of the results presented in this 
paper. The case of quantum dot with a non-collinear mag- 
netic order has been discussed in fl2l . Before addressing 
the case of magnetic rings, we discuss several fundamental 
points. In addition to the effects of electrons on the DWs 
we are also interested in how the ferromagnetic nature of 
the system and the presence of the DWs affect such well- 
known phenomena as the Aharonov-Bohm effect 11311141 . 
persistent currents U5U16I1 . and conductance fTTll . 

It is established that strong carrier scattering and in- 
terference results in long-range interactions between im- 
purities on metal surfaces. The question of whether and 
how the carriers' spin dependent scattering mediates inter- 
actions between DWs is addressed here. We have found 
that the DWs are coupled thusly: due to the scattering from 



1 This can be explained by simple magnetostatics accounting 
for the demagnetization fields. 
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the first DW a spiral spin density wave is created. This then 
acts as a spatially non-uniform torque on the second DW 
whose energetically stable shape and position therefore has 
a non-uniform dependence on the distance from the first 
DW. This is different from the ordinary spin-torque trans- 
fer in bulk spin valve systems IU8lll9ll20l or magnetic tun- 
nel junctions 12111221 insofar as in this case the DWs spatial 
arrangement, in addition to the magnetization direction, is 
current controlled. 

There are a number of experimental works on permal- 
loy ferromagnetic rings (8]|9l. In these works the exis- 
tence of the Aharonov-Bohm effect was experimentally 
confirmed 12311241 and the studies of the magnetoresistance 
effect have also been presented 1125112611271 . However, it is 
still not clear how the variation of magnetization in these 
systems affects the state and dynamics of the DWs. Though 
our aim is to consider DWs in the nanoscale rings, we dis- 
cuss first the rings without DWs. The main question is how 
the vortex state in a magnetic ring affects its electric prop- 
erties. In particular, how the corresponding edge states look 
like and how persistent currents and Aharonov-Bohm type 
effects are modified by choosing different magnetization 
configuration. 

The first question to answer involves analyzing the 
electron wavefunctions of the ferromagnetic ring II28II291 . 
Note that in the following we call the magnetization of 
ring "ferromagnetic" when the magnetization vector is di- 
rected along the 9 angular unit vector in the cylindric coor- 
dinates of the ring. In this configuration there is no energy 
cost from the demagnetization field, unlike the true ferro- 
magnetic case. The vortex state instead incurs the smaller 
penalty of not having exactly parallel magnetic moments 
everywhere, though locally they will be approximately par- 
allel for a large ring. Note also that as we are dealing with 
a ring, there is no energy cost due to the vortex core. 

The Zeeman splitting of the electron spectrum in a 
magnetic field has been included in Refs. 130113 1 1 . The 
overall effect of this, as far as the Aharonov-Bohm effect 
is concerned, is to alter the dephasing length scales. The 
spin-orbit coupling has also been included but only per- 
turbatively where it mixes the spin channels. Though the 
spin-orbit coupling could be considered analogous to what 
we wish to calculate (which also couples angular orbital 
momentum and spin) we solve this model exactly. 

Another interesting idea concerns the magnetically tex- 
tured mesoscopic ring of Ref. 113211 . In this model, one con- 
siders the inhomogeneous magnetic field distorted from the 
z-direction around the ring, like a crown. This induces a 
form of spin-orbit coupling which can lead to persistent 
currents even in the absence of a magnetic flux threading 
the ring. The method can be useful for decoupling the or- 
bital and spin degrees of freedom by using a Feynman path 
integral representation. As far as the domain wall motion 
is concerned some efforts have been also made for semi- 
rings l33l . It has been also proposed to take into account 
the inhomogeneous anisotropy, which arises naturally in 



magnetic nanoring [ 34 . 35 . 35 1 . It was shown that the effect 
of the resulting Berry phase of electrons is quite similar to 
that of an external field affecting the equilibrium magneti- 
zation profile. 

2 Wires 

2.1 Theoretical model We consider first a long mag- 
netic wire with two DWs, down which a current, /, is 
passed (as shown schematically in Fig. [TJ. Assuming that 
the distance zq between DWs is larger than the phase co- 
herence length Lj„ we can consider them as two indepen- 
dent scatterers, but for zq < the current transmission 
mediates DW coupling. For definiteness, we assume that 
one of the DWs (located at z = 0) is pinned, e.g. by a ge- 
ometric constriction, and concentrate on the effect of the 
current on the second DW, initially (i.e., for 1 = 0) lo- 
cated at z = zq. For 1 = each DW has an extension 
L, We model the one-dimensional wire with the following 
Hamiltonian, H, of noninteracting electrons coupled (with 
a coupling constant J) to a spatially non-uniform magneti- 
zation (DW) profile M(z): 



H 



dz a) a (z) 



2m 



jcr af) -m{z) 



ap{z), (1) 



(we use units with U = 1). Here a^ a and a a are the cre- 
ation and annihilation operators of electrons with spin a. 
We use a local unitary transformation T{z) 136,37,38] 
(often referred to as a gauge transformation) which trans- 
forms the second term with nonuniform magnetization into 
the constant Zeeman splitting term and an additional spin- 
dependent spatially varying potential U a p(z). For kpL > 
1 this potential can be treated perturbatively [36 37,38,391 
(k,F is the Fermi wave vector) whereas the sharp domain 
walls, i.e. when kpL < 1, require a different formalism 
HQ). 

The matrix T(z) is defined by T^(z) er ■ n(z) T(z) = 
a z , where n is the unit vector along M, M(z) = M n(z). 
Then the transformed Hamiltonian H = T^(z) HT{z) is 



H 



dz a' a (z) 



2m 



+ U aP (z)- JMa z ap 



ap{z). 



(2) 



For a wire with two DWs we can parameterize the magne- 
tization profile by two angles <p(z) and 9(z) (cf. Fig.[TJ 



n(z) = (cos# sin ip, sin (9 sin ip, cosip), 
z 



(3) 



ip(z) 



tanh [— ] ) + cos" 
L 



i( tanh (4) 



=—tp x (z) =— ¥>a(*) 

(for details see Ref. ITOl ). The relative orientation between 
the two DWs (situated, respectively, at z = and z = zq) 
is set by angle 9(z). We set 9\ = at the first DW and 
#2 = 8q around the second one (see Fig. Q]). In the fol- 
lowing we consider the case of zq > L, which can be 
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treated perturbatively as the coupling of the DWs is rela- 
tively small. Then we can write U(z) w U\(z) + ^2(2), 
where ( j = 1 , 2) 



U,(z) 



8m 



<(f) 
4m 



<p'!(z) | ^{z)d, 
4m 



2m 



cos 9 j(z) 



2m 



1 



(5) 



Note that this approach is generalizable to any number of 
DWs which are sufficiently far apart. 

As shown in Refs. II 3 7 II 3 8 11391 , for a single DW and 
kpL > 1 a perturbative approach is appropriate for treat- 
ing the electron scattering from the DW potential (01. Thus, 
assuming ip°i z ) to be the wave function of electron with 
energy e in the wire without potential (0, we find the first- 
order correction due to the perturbation U\{z), i.e. due to 
scattering from the first DW, as 



dz' G e (z,z')C/i(z')V' (2')- 



(6) 



The Green's function G e (z, z') corresponds to the unper- 
turbed Hamiltonian (2) with U(z) = 0. It is diagonal in the 
spin space with elements 



G £(T (z, z ) — — 



where k„ 



k° + t^S- for lifetimes r a > e. 



(7) 



and 



fc£ 4 = [2m(e + n± JM)} 1/2 . 

The interaction energy of the two DWs due to the single 
scattered state tp su (z) = ^ a {z) + 5ij) Ea {z) is 



AE„ = 



dz 8\l)l a (z) U 2 (z) 5i)) ea (z). 



(8) 



To find the total interaction energy we should sum up the 
contributions of all scattering states in the energy range 
between ef and £p + eAip/2, for an applied voltage 
eAtp/2 <C ep. Then we find the current-induced coupling 
of the DWs: 



AE 



eA(j> ( AEf 

72V v 



AE } 



(9) 



where v a = k®/m is the velocity of spin up and down 
electrons at the Fermi level. 

For numerical calculations we use the parameters of 
magnetic semiconductors IfT PD . It should be noted that in 
the case of metallic nanowires the ID limit is also achiev- 
able l42l . Here we take the parameters as in Ref. 0, i.e. 
Xp = 6 nm; an effective mass of m = 0.5m e (m e is free 
electron mass); L = Xp; JM = 15 meV; £p = 83.7 meV; 
and eA(j> = O-le^. We also assume relatively large mean 
free path of I = 500 nm, corresponding to essential order- 
ing of Mn ions. The width of the wall may be as small as 
the atomic size in the presence of constrictions 143 441151. 
hence well below the DW lengths in bulk materials. In such 
a situation, the interaction between the DWs increases due 




AE(zu,9o) 0.00 
meV 

-0.05 



Figure 1 (Color online) Top panel: A schematic showing 
the DWs magnetization profile (thick arrows). L is the DW 
width, zq and 6q are respectively the DW position and ori- 
entation with respect to the DW at z = 0, J is the current 
direction. Lower panel: Interaction energy AE(zq, 9q) as a 
function of zq and #o- The solid curve is for zq = 30 nm, 
the dashed curve is for zq = 37.5 nm, and the dotted curve 
is for zq = 45 nm. 



to strongly enhanced DW scattering I45ll46ll47ll48l . The in- 
teraction energy shown in Fig. [TJ depends periodically on 
the DWs relative angle 9q and distance zq. It results in an 
oscillating motion of the DW along the axis z as well as 
an oscillating direction of DW polarization. In its turn this 
will have an effect on the spin torque driving the DW dy- 
namics, which we will now focus on. 

We find the spin density due to a single transmitted 
electron wave of spin a is 



S a (z) = iPl(z)T(z)<TT\z)TP Ea (z), 
and the total current-induced spin density is 



(10) 



S(*) = f^ + ^). (ID 

We find that the correction to the spin density follows the 
magnetization profile with additional Friedel oscillations, 
which are a superposition of two waves with periods kpi 

and kp, . The oscillations in the spin density are smaller in 
magnitude than the overall spin density profile and decay 
with increasing z. 

We calculate the current-induced torque acting on the 
magnetization at z, for the second DW located at zq, from 



AT{z 7 z ,9 ) 



-M(z,z ,e ) x AS{z,z ,9 )Al2) 



where 7 = <?^b, g is the Lande factor and /is is the Bohr 
magneton. We assume a thin nanowire with a cross section 
of a cs = 100 x 20 nm 2 as in Ref. Q. In Eq. (15) AS 
is the correction to the electron spin density due to scat- 
tering. The calculated torque on the second DW is shown 



Copyright line will be provided by the publisher 



4 



Nicholas Sedlmayr et al.: Non-collinear magnetization 




Z-Zq 



Figure 2 (Color online) The x-component of the current- 
induced spin torque, as defined in Eq.[l2] acting at the sec- 
ond domain wall as a function of z and 8q. 
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Figure 3 (Color online) The magnetization dynamics of 
the second DW when the anisotropy is included, see 
Eq. [13] The solid curve is the x-component, dashed the 
y-component, and dotted the z-component. Taken at the 
centre of the domain wall. 



in Fig. |2] for the x-component, where zo = 50 L and 
M « 5.56 x 10 4 Am _1 were used BTI . The correction to 
the spin torque shows that the force upon the DW depends 
strongly on their relative orientations. 

To inspect the current-induced dynamics of the DW at 
z = zq, we evaluate the accumulated spin density that acts 
on the DW at z = Zq. The DW magnetization dynamics 
are then modeled using the Landau-Lifshitz equation with 
the effect of magnetic anisotropy and damping: 

d t M = — — M x S[M1 + ^ArM x xM x 

a cs 

+ -^M x <9 t M. (13) 

As an initial condition we assume that the magnetization 
profile in the wire without electric current is described by 
Eq. We should note that the relative orientation of the 
walls at the start of motion does play a role in the type of 
motion we see. Here we present it for an arbitrary configu- 
ration. 

We take the anisotropy constant to be K' = — 10, and 
therefore the x-axis to be a hard magnetization axis. Figure 
[3] shows the effects of anisotropy on the domain wall mo- 
tion. The anisotropy dampens motion in the x-direction, 
thus exacerbating the y and z oscillations. This is also in 
contrast to the case where we ignore the first domain wall. 
In this case, although the anisotropy does introduce motion 
around the centre of the domain wall it does not involve a 
decaying ^-component, see Fig. [4] 

Finally let us include weak Gilbert damping also (a — 
0.01), we then find a switch over between behaviour dom- 
inated by the anisotropy and behaviour dominated by the 
damping. See figure|5] 

3 Quantum rings 




0. 5.xl0~ 7 1.x 10~ 6 

*(s) 

Figure 4 (Color online) The magnetization dynamics of 
the second DW without the influence of the first DW, when 
the anisotropy is included, see Eq. [13] The solid curve is 
the x-component, dashed the y-component, and dotted the 
z-component. Taken at the centre of the domain wall. 




Figure 5 (Color online) The magnetization dynamics of 
the second DW when both anisotropy and damping are in- 
cluded, see Eq. [13] The solid curve is the .x-component, 
dashed the y-component, and dotted the z-component. 
Taken at the centre of the domain wall. 
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3.1 Model The following Hamiltonian 

2 



H 



1 

2 m 



-drirdr) +\(d e - i— 
r r z \ 0o 

-J Meg ■ a + V(r), 



(14) 



describes a two dimensional ferromagnetic ring of elec- 
trons confined by some radial potential V(r). Here we use 
the usual notations: J is the exchange coupling constant 
and m is the electron mass. In correspondance to the ring 
geometry it is convenient to use the polar coordinates. We 
denote 



= d> A • dr = 

J 1 



B • dS, 



(15) 



(fro = h/eis the magnetic flux quantum, B is the mag- 
netic field and A is the vector potential. One can simplify 
Hamiltonian (14) making a local rotation in spin space to 
remove the 9 dependence from the JM term. We use the 
following unitary transformation: 



1 



(16) 



Note that this transformation preserves the same boundary 
conditions for the wavefunctions. The transformed Hamil- 
tonian is H' = TJT (6) HXJ (9) and the corresponding wave 
function is ifr(r, 8) = U'(6)ij)(r, 0). We obtain 



H' = - 



1 

2m 



d r (rd r 



1 



do 



■4^2- ^1 d o- 1 — 



where we introduced 



- J Ma, 



.<fr' 



V(r) 



= 4> — <fro/ 2. The latter means 
the appearance of an additional phase of an electron </>o/2 
related to the ring geometry. 

The angular dependence can be dealt with by a Fourier 
transform. We use the following definition: 



V>(r,0) = e m if>i(r) and 



l=- 

,■2 TV 



1 rlf) 

Mr) = J o ^e- iW ^(r,6). 



(18) 



Finally we obtain the following transformed Hamiltonian: 



2m 



1 



1 

-J Ma, 



I', 



+ V(r), (19) 



where l'^ = I — (fr' /(fro- [From this point one can find the ap- 
propriate one-dimensional Hamiltonian |50| which can be 
further used to calculate the spin current in the ring without 
the DWs. 




Figure 6 The persistent current, j e {4>), in a one- 
dimensional simple ferromagnetic ring. 



3.2 Onions and vortices Once we understand the 
problem of the vortex state without the DWs, we can con- 
sider peeling back the layers of the onion state. The onion 
state refers to a ring in which there are two domain walls 
present, and hence two regions of oppositely directed mag- 
netization. One can consider such state in the limits of 
sharp and adiabatic DWs (as compared to the electrons' 
wavelength). The following refers to the ID model. Note 
that one should be careful in taking the ID limit when a 
spin-orbit interaction is present [1501. 

3.3 Magnetic ring without the domain walls In 
this case we can make the transformation of Eq. ( fT6l >. Then 
using H' = W(9)HV(6) and tp(6) = Ut(0)<0(0) we 
obtain 



(17) H' 



1 



2mp 2 
1 

v 2mp 2 



oe - 1 — 



1 1 



do - t— 



4 2mp 2 



JMa z 



(20) 



After the Fourier transform we diagonalize the Hamilto- 
nian; H" = W(l)H'N(l) andt/>(Z) = Nt(Z)^(Z). 



N(Z) = 



1 



1 iA2 



Af W 



(21) 



Ay = -il'^/iJv^ + p 2 - m) with fi = 2JMmp 2 . Then 
we obtain the exact formulas for the energy spectrum 
1 

4 



£/(i,2) 



"F I A* 



2 m p 2 
2 



1 + A+ 

\-Ai 



2ilL 



1-A 



(22) 



Now we can directly calculate the persistent current. 
Using 



(23) 
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Figure 7 A schematic view of a ferromagnetic ring with 
an applied magnetic field. 



(where i = 1, 2) we define the charge current as j e {<j>) = 
ji(4>) + j2{4>)- It is presented in Fig.|6]as a function of <f>, 
3.4 The onion state with adiabatic domain walls 

In this section we use a modified version of our model for a 
wire with two domain walls ifTUl shown in Sec. 2. By map- 
ping the circle onto a line, and directing the magnetization 
along eg instead of e z , we can easily construct a suitable 
model. The mapping is given by 



x(9) = y tan 

lJ 



(24) 



This mapping distorts the relative sizes of the domain walls 
depending on where on the ring they reside, and its accu- 
racy increases with increasing the ring diameter. 



cp(9) 



- 1 [tanh[a:(0) - x(0i)]] 



7T — COS 



=<pi(e) 

[ta,nh[x(6) - x(6 2 )]] 



(25) 



= <P2(6) 



where 0j is the position in the ring of the zth wall, which 
we take to be at an angle rji to the z-axis, see Fig. [7] The 
magnetization is then 



M(0) 
M 



= cos[tp(9)]eg + sin[?/(0)] sin[</?(0)]e r 



+ cos[r/(0)] sin[yj(0)]e z , 
and the Hamiltonian is given by 
1 



H = 



2mp 2 



- JM(0).<7. 



(26) 



(27) 



where 7/(0) is a function, which takes the value 7/j around 
the zth wall and slowly interpolates between the walls. 



We can locally rotate the spin direction as before. The 
resultant Hamiltonian, 



H' 



1 



2mp 2 



d e - i- 



1 



i [de - i- 



4 \ (po, 

- JM(0) • cr, (28) 
includes the following magnetization profile: 



M(0) 
M 



cos[^(0)]e 2 



- sin[7/(0)] sm[cp(9)]e x 
-cos[r/(0)] sm[ip(9)]e y . 



(29) 



After that we can make a gauge transformation to reduce 
the magnetization profile to a constant Zeeman term and an 
additional potential. However, the non-diagonal portions of 
the transformed Hamiltonian complicate significantly the 
exact calculation of the energy spectrum. 

Therefore, we consider a simplified approach which is 
justified for large magnetic ring, kp ^> 1, where k is the 
electron wave vector along the ring. Using Eq. (29) we 
write 



H' = --f- 



1 

2m 



1 



l -\d z 



Kp 



- JM(0) ■ cr, (30) 



where z = p9 is the coordinate along the wire. Due to 
the periodicity in angle 9, the condition of periodicity z —> 
z + 2irp should be preserved. If kp 3> 1 the second term 
in Eq. (30) is small and can be taken into account perturba- 
tively. If we neglect the small correction we come to the ID 
Hamiltonian, which can be treated using the same methods 
as the linear model of Sec. 2. The main difference is that 
the gauge potential related to the DWs is an additional cor- 
rection to the gauge related to the ring curvature and to the 
magnetic flux through the ring. 

It should be noted that the ring topology allows only 
an even number of DWs. In the case of two DWs on the 
ring, the current-induced interaction between them can be 
described by the spin torque formula Eq. (fT2l . On the other 
hand, the persistent current in the magnetic ring, Eq. ( 1231 , 
would necessarily induce the current-induced force and an 
interaction between the walls, which puts them both into a 
kind of circular motion. However, one can also expect that 
the damping associated with the DW motion would lead 
to a slow decay of the persistent current. We can assume 
that a slight pumping of power from an external source can 
make the motion of DWs on the ring non-decaying, which 
makes this problem attractive for possible applications. 

4 Summary We have shown how the presence of 
more than one domain wall in a wire affects the domain 
wall dynamics and the energy of the electron system. We 
calculated the energy due to the current mediated coupling 
between the domain walls and demonstrated how the do- 
main wall dynamics is modified by the current passing 
through a pair of the domain walls in the wire. 
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Then we consider ferromagnetic nanorings with two 
domain walls. In this case the DW dynamics is also af- 
fected by the coupling between the walls. Besides this, 
there appear additional contributions to the gauge poten- 
tial related to the ring geometry of the ferromagnetic state 
and to the magnetic flux through the ring. In the case of 
a relatively large ring with two domain walls, the problem 
can be solved as for the wire with the periodic condition 
for eigenfunctions, resulting in the quantization of values 
of the wave vector k along the ring. 
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